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Abstract 

We study the mass, width and couphngs of the hghtest resonance multiplet with 
/ (J^^) = 1 (1 ) quantum numbers. Effective field theories based on chiral symmetry 
are employed in order to describe the form factor associated with the two-pseudoscalar 
matrix element of the QCD vector current. The bare poles of the intermediate reso- 
nances are regularized through a Dyson-Schwinger-like summation. We explore the role 
of the resonance width in physical observables and make a coupled-channel analysis of 
final-state interactions. This provides many interesting properties, as the pole mass 
MP°'° = 764.1 ± 2.7 tti MeV. At energies E >1 GeV, a second 1(1"") resonance 
multiplet is considered in order to describe the data in a more consistent way. From 
the phenomenologically extracted resonance couplings, we obtain the chiral coupling 
Llifio) = (7.04 ± 0.05 • 10"3, at fio = 770 MeV, and show how the running with 
the scale fi affects the resonance parameters. A I/Nq counting is adopted in this work 
and the consistency of the large- Ac expansion is tested. Finally, we make an estimation 
of the contribution from diagrams with resonances in crossed channels. 



1 Introduction 



It has become evident that Quantum Chromodynamics (QCD) is the correct theory which 
describes hadronic processes [0. In the high-energy region (i? ^ 1 GeV) the theory accepts 
a perturbative description and, accordingly, many calculations up to several orders in the 
perturbative expansion parameter have been performed. These theoretical results have 
been successfully tested in many high-energy experiments. Nevertheless, since the running 
coupling constant as(yu) increases as the energy decreases, the perturbative expansion in 
powers of Us breaks down at energies E 1 GeV. In this paper the problem of describing 
the E <1 GeV region by employing effective theories of QCD [0, ^ will be analyzed. 

When we study processes at energies much lower than the heavy quark masses, the degrees 
of freedom corresponding to heavy quarks decouple ^ and QCD, with only the light quark 
fields, yields a proper description. In the massless limit, the QCD lagrangian shows chiral 
symmetry: the left-handed and right-handed quark fields can be rotated independently under 
the SU{nf)L ® SUinf)^ flavour chiral group, where nj is the number of light quarks. The 
symmetry is spontaneously broken to the SU{nf)y subgroup and nj- — 1 massless Nambu- 
Goldstone bosons appear, associated with the broken generators. Nonetheless, as the light 
quark QCD lagrangian has small non-zero mass terms, chiral symmetry is also broken explic- 
itly and the Nambu-Goldstone bosons gain small masses. These bosons have = 0~ and 
are identified with the triplet of pions, in the SU{2) case, and the (vr, rjg) octet of light 
pseudoscalars for SU{3). 

The low-energy chiral effective field theory describing the dynamics of the lightest pseu- 
doscalar multiplet was first developed for the SU{2)l® SU{2)r symmetry group 0, and was 
later generalized to the three flavour SU {?>) l® SU {?>) r case [|]. We will use the latter in order 
to include kaon interactions in our study. Chiral Perturbation Theory (xPT) P, H, 0, ||, 
describes the physical low-energy amplitudes as an expansion in powers of quark masses and 
momenta over a characteristic chiral scale A-^ ~ ^T^f-K ~ 1 GeV, with f^^ = 92.4 MeV the pion 
decay constant. 

The expansion in powers of momenta over A,,, deteriorates as the energy of the process is 
increased and, in order to reach the relevant accuracy, one needs to add higher and higher 
chiral orders to the xPT lagrangian. In the resonance region one must introduce a different 
effective field theory with explicit massive fields to describe the degrees of freedom associated 
with the mesonic resonances. In the eighties, Gasser and Leutwyler worked out an SU{2)l ® 
SU{2)r lagrangian describing the pions and the vector resonance p(770) Later on this 
work was extended to the nf = 3 case |jlO|, developing the Resonance Chiral Theory (R^T). 



Further studies on the R^T and xPT lagrangians constrained the resonance chiral couplings, 
employing the QCD short-distance behaviour of appropriate Green functions . 



Once the resonance fields are explicitly included in the effective lagrangian, the chiral 
counting becomes ineffective because the masses of these resonances are of the same order 



1 



than the chiral characteristic scale A^,.. However, an expansion of QCD and its low-energy 
effective field theory in powers of I/Nq, with Nc the number of quark colours, appears to be 
suitable |]12|. In the large- A'",^ limit, the hadronic description reduces to tree-level processes 
without hadron loops. As the 1/Nc expansion seems to yield a proper description of Nc = 3 
QCD, it seems also appropriate to expand the R^T results in powers of l/Nc [|13|- To a 
certain extent, this reduces to just counting the number of loops. 

At leading order (LO) in 1/Nc, RxT yields a good description of many phenomena. 
However it fails when the energy approaches the bare mass of a resonance. This situation is 
common to every unstable propagating state in a Quantum Field Theory when its propagator 
turns on-shell. It is solved by the Dyson-Schwinger summation of one particle insertion 
blocks (IPI), which provides the unstable particles with an imaginary absorptive part in the 
resonance propagator. This summation must also be done in R^T, with some prescriptions. 



but essentially in the same way. In Refs. |]T4|, the p-channel was studied and an appropriate 



off-shell width for the p(770) resonance was obtained. 

In this paper we continue the work put forward in Ref. [|15], extending it to a coupled- 



channel analysis. We will study the vector form factor (VFF) jn, |T^, [TB|, |T^, |TB| and overview 
the correlator of two QCD vector currents and the corresponding partial-wave scattering 
amplitude. It will be shown that our coupled-channel description of the resonance width 



agrees with the one in Ref. [|I5|, obtained with a single-channel treatment. From the Dyson- 
Schwinger summation, we find that the rescattering dresses the bare propagator in a universal 
way. The induced correction only depends on the intermediate IPI blocks and not on the 
final or initial states of the process. 

We briefly describe the basic ingredients of the R^T effective action in Section 2. The 
Dyson-Schwinger analysis of the different observables is performed in Section 3. In Section 4 
the obtained results are matched to the C(p^) xPT description and compared with the data 
in Section 5. The xPT coupling Ivg(/i) is calculated here and a test of the l/Nc expansion is 
also performed. The small corrections induced by resonance exchanges in the t-channel are 
estimated in Section 6. Our conclusions are finally given in Section 7. Some technical details 
have been relegated to the Appendices. In particular, a generalized formal summation of 
diagrams with two-body topologies is presented in Appendix C. 



2 Resonance Chiral Theory 



We will work with the SU (3) octet of light pseudoscalar bosons, interacting through the 
;^PT lagrangian |, 0: 

4?T = y + (1) 
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where (...) is short for the trace over flavour matrices and / ~ /yr is the pion decay constant 
at lowest order. The tensors = i[u^ {d'^ — ir'^)u — u{d'^ — il^)u^ and x+ = tt^X"^^ + "^X^"^ s-'^^ 
functions of the left, right and scalar external sources Z^, and x 0]- The pseudoscalar 
fields 

„0 



m 

K- 



I '18 

'V2_J v/6 



2m 



(2) 



are parameterized through the SU{3) matrix u = exp 



fV2 



£28 / 



The interactions of the Nambu-Goldstone bosons with the lightest multiplet of vector 
resonances 



V, 



V2 



K* 



K 



*o 



(3) 



are given 



by iig 



RxT 



(4) 



where = u F^^u^ + u^F'^^u, with -F^'^ the field strength tensors of the left and right 
external fields ||5|, |^. We use the R^T lagrangian in the antisymmetric formalism provided 
in Ref. It was demonstrated in Ref. |11] that this antisymmetric description of the vector 
fields is equivalent to the more usual Proca formalism plus the 0{p^) xPT lagrangian with 
its couplings Li constrained by the short- distance QCD behaviour. 

In general, one should consider a set of vector resonance multiplets VjiJ with couplings 
Fy^ and Gy.. At low energies (a/s < 1.2 GeV), the lightest multiplet yields the dominant 
contributions. However the tail of the second nonet may generate sizable corrections which 
must also be taken into account in the ^/i ~ 1 GeV region. 



3 The Vector Form Factor 



Let us consider the hadronic matrix element corresponding to the production of two pseu- 
doscalars with I = J = 1 through the charged d'y^u vector current: 

( P-iPi) P\P2) I d^'u |0) = V2{p,- p^r :F^''\q') , (5) 

with q = pi + P2- The label P denotes the pair of pseudoscalars which are produced in the 
final state, either vr^vr^ or K~K^ . The Lorentz structure is fixed by current conservation in 
the isospin limit. 
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Figure 1: VFF at leading order in 1/Nc with one and two vector resonances. For the two-resonance 
case we have adopted the input parameters Mvi = 775 MeV, Myj = 1450 MeV, Fy-^Gvi//^ = 1.1 
and Fy^GyJf = -0.1. Data from ALEPH p| and NA7 pO]. 



At leading order in 1/Nc, the vector form factor TQ^\q^) is easily computed through the 
diagrams shown in Fig. ||(a). We put together the two J^lf\q'^) functions in the vector 

The requirement that the vector form factor should vanish at infinite momentum transfer 
constrains the resonance couplings at LO in I/Nq to satisfy the short- distance QCD rela- 




tion 11, 13 



i-E^ = o. (7) 

i J 

If only one vector multiplet is considered, then Fy^Gy^j p = 1 and one gets the familiar 
vector-meson dominance expression 

M^') = f _i ) • (8) 

My^-q V V2 J 

This yields a rather good description of the data in the region < 0.7 GeV, below the 
p(770) peak. Chiral loop corrections are subleading in the 1/Nc counting and turn out to 
be rather small in this case. Other resonances can also be included. The relevance of the 
large-iV(7 expansion to approximate the physical vector form factor is clearly seen in Fig. |l|, 
either with just one resonance or including a second multiplet. 

The vector couplings are as well constrained in the large-A^^ limit by the relation 

E2 Fy. Gy. — F'yr, , . 
^-^^ ^ = 0, (9) 

I Vi 



4 




(b) 

Figure 2: Effective vertices for the vector current insertion producing two pseudoscalars (a) and 
for the two pseudoscalar scattering (b). The first terms come from >C^pt and the second ones from 
the interaction via an intermediate resonance due to the CYrj, lagrangian. 



provided by the short- distance QCD conditions over the axial form factor |jTT], [13 



For the simplest situation with a single resonance exchange 

V2 / in the large-iVc limit |ll| 



constraints yield Fy^ = 2G 



Vi 



the short- distance QCD 
However, since we are going 
to work at higher orders in 1/Nc, we will leave these couplings free and will test afterwards 
the deviation of their experimental values from the large-iVc predictions, that we expect to 
be small. 



3.1 Dyson- Schwinger Summation 

At energies close to the mass of a resonance we need to know the denominator of the 
resonance propagator beyond the leading, bare, order in l/Nc- What is usually done is 
a Dyson-Schwinger summation, as for instance in the QED photon polarization. That is, 
summing diagrams composed by a series of propagator, IPI block, propagator, . . . , and so 
on. This summation regularizes the pole of the bare propagator. It gives a self-energy with 
its corresponding absorptive part, up to the perturbative order employed for the IPI block. 
In RxT, however, at the same order than the resonance-exchange contribution there is also a 
local interaction from the C^^pt lagrangian. The Dyson-Schwinger summation must be then 
slightly modified. One constructs effective current vertices and effective scattering vertices 
15| , by adding the contribution from intermediate resonance exchanges in the s-channel 



(2) 

to the local xPT interaction C^Jt. Both contributions are of the same order in the 1/Nc 
counting. These effective vertices, shown in Fig. are independent of the explicit formulation 



adopted for the spin-1 fields [|TT], |T5|]. If we use the Proca formulation we have to take into 
account the local interaction from the 0{p'^) xPT lagrangian as it is described in Ref. |jlT 



The inclusion of the local vertices is not important on the resonance peak but it turns to be 
relevant away from it. 
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Pi' P2' P,' Pn' 

Figure 3: Diagrammatic summation at N loops. 



For the moment, we are only interested in the imaginary part of the self-energy. Therefore, 
we will concentrate in the sum over diagrams with absorptive cuts. For the range of energies 
we are interested, the most relevant contributions come from intermediate states with two 
pseudoscalars; states with a higher number of particles being suppressed by phase space and 
chiral counting. Thus, we are going to sum diagrams^ constructed with an initial effective 
current insertion connected to an effective scattering vertex through a two-pseudoscalar loop. 
The pair of outgoing pseudoscalars from the scattering vertex are again connected to another 
effective scattering vertex through another two-pseudoscalar loop, and so on, as it can be 
seen in Fig. |[ 

The off-shell effective current vertex shows the momentum structure 

v,^ = V2[foPr + f'oPr] (Pl-P2)^ m 

with Pj^'^ = g^^ — 2^ and = the usual transverse and longitudinal Lorentz projec- 
tors. In the isospin limit, the second term with JF'q vanishes when the outgoing pseudoscalars 
are both on-shell. Notice that this off-shell function depends on the adopted parameterization 
of the fields, but the final on-shell amplitude does not depend on it. 

When the current insertion Vq is connected to a successive number of loops and effective 
scattering vertices one gets = TmPt" + ^' NPf (Pi ~ P2Y , where N is the number 
of intermediate loops in the diagrammatic chain shown in Fig. ^. Thus the momentum 
structure remains. Inductively, from to iV -|- 1 loops we can observe the linear recurrence 
T^^^i = J2j ■M^-' where i = 1 stands for tttt and i = 2 for KK. This feature can be 
expressed in the matrix form 

fN+i = MfN = = ■■■ = M^'^'fo, (11) 

with jF/v the vector form factor at loops. The 2x2 recurrence matrix takes the form 

M = -^-'Tl^^-' {192nB22), (12) 

^ This diagrammatic construction solves the Bethe-Salpeter equation in an iterative way. The effective 
vertices provide the corresponding "potentials" at LO in 1/Nc- 
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with the diagonal matrix S = diag (cr^, cr_ft-), being crp = ^1 — Aim^p/q^. The matrix 

T^^5iF{^-?^MS^^- 1 . rM^. (13) 

is the s-channel partial- wave scattering amplitude with / = J = 1, at LO in 1/Nc [Fig. 0(b)]. 
The diagrams with resonances in the crossed channels produce a tiny contribution which 
will be taken into account in Section 6. We can also observe in Eq. ( |T^ the diagonal 
matrix = diag 522'* ; B^^^, with the two-propagator Feynman integral b'^^ given in 
Appendix A. 

Summing the result in Eq. (|TI|) for any number of loops, one gets a geometrical series 
which can be easily handled: 

OO / OO \ -1 

^ = E -^^ = E -^"^ -^0 = {l-M)-' = ^— TTTY^o. (14) 



Af=0 



l-ii{M) 



The last identity is not trivial. The S ^T^qS ^ matrix is proportional to a dimension-one 
projector and jFg is eigenvector of this projector. Thus, acting over JFq reproduces again 
the vector times a number. The mathematical details can be found in Appendix B. 

Afterwards a more complete calculation of the form factor will be developed. At the 
moment only the absorptive diagrams have been included and only the imaginary part is 
under control. Moreover, let us consider the simplest case of a single resonance exchange. The 
factor together with the initial jFg generates a complex denominator My^ ~ (f — ^{(f')'- 

a non-controlled real part plus a well defined imaginary term, given by 

Ime(g2) = Im^^(g2)+Imei^(g2) . (15) 

The bubble loop summation provides an imaginary contribution which gets separate con- 
tributions from the tttt and KK channels. The corresponding partial widths are provided 
by 



Im ip{q') = Cp (^M4 + ^ Q') Oiq' - 4m^) , (16) 

with = 1 and Ck = \- When we substitute the coupling at LO in l/Nc, Gvi = f/V^, 
these imaginary terms Im ^p{q'^) agree with the partial widths MpT^p^\q'^) obtained in 



Refs. [14, 15] from a simplified single-channel analysis. This energy dependence for the 



width was long ago considered by Gounaris and Sakurai from general arguments [^. As 



well, they had exactly the same logarithm in their work that the one which naturally appears 

(p) 

in our calculation of the absorptive contribution through the Feynman integral i?22 • 

The correlator of two vector currents and the I = J = 1 partial-wave scattering amplitude, 
can be computed in a similar way. For the correlator we begin with a current effective 
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vertex, like for the form factor, and connect it to a A^-loop final-state interaction, which 
ends into another current effective vertex. For the scattering amplitude we start from a 
scattering effective vertex and go on connecting loops and scattering vertices in the same 
way. A similar rescattering effect appears in the three quantities. The resulting (s-channel) 
I = J = 1 scattering amplitude takes the form: 

^ = S ( E -^'') = S (1 - M)-' Tl, = Y—^lM} ■ (17) 

The matrix structure (1 — only depends on the scattering effective vertex and on 

the two-intermediate particle loop. As these are identical for the three quantities (VFF, 
correlator and scattering), the final-state interaction dresses the bare resonance pole in a 
universal way, providing the same complex pole for all processes. 



4 Low-Energy Matching Conditions 



All the former calculations must reproduce the QCD low-energy behaviour provided by the 
xPT framework. This allows to fix the polynomial ambiguities at a given order in the chiral 
expansion. We can identify the momentum expansion up to 0{E^) of the resummed vector 
form factor (|l4|) with the standard 0{E'^) xPT calculation in the usual MS — 1 scheme [23|. 



At leading order in l/iVc, we have the well-known relation |T0|, [13 



^ 2M? ' 



Vi 



with Fy-, and My- the bare parameters of the R^T lagrangian. 



Keeping 1/Nc corrections, the 0{E ) matching determines the regularized function B22 



r,{P) 



up to the considered chiral order, to be uM, 15, 1 



1 



1927r' 



cr". 




5 
3 



> rrip 



(19) 



where SLl{fj,) = L^fx) — Lg\j^^_^^. The renormalization scale dependence of the O^E"^) 
XPT coupling Lg{fi) cancels out with the term In {mj,/ fi"^). The resulting vector form factor 
from (|l4f) takes then the form 



1 + E 



FyG 



Ml 



1+ 1 + E 



\ V 

P Ml-q^ P 



B. 



r,{iT) 
22 



+ 2^22 



r,{K) 



1 
V2 



(20) 
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With the information obtained from the VFF we obtain as well the (s-channel) I = J = 1 
partial-wave scattering amplitude, 



T 



96nf' 



( 



a: 



1+ i + E 



q^ 



2q' 



B. 



r,{iT) 
22 



+ 2^22 



r,iK) 



■ V2 



' V2 



a 



, (21) 



with -822^^ being the same than in the VFF due to the optical theorem. 



4.1 Scale Running 



When the low-energy matching was performed, the unfixed (5Lg(/x) parameter was left. It 
appeared as an extra constant in 822^^ ■ This also pointed out an ambiguity in the election 
of the scale and in the renormalization scheme, usually MS — 1 but not the unique one. For 
simplicity we will analyze this feature in the single resonance case and with the leading values 
of the couplings, Fy^ = V2f = 2Gvi- this situation the VFF, for instance, becomes 



P 



^22 



2 



r,iK) 
22 



1 
V2 



where 



P 



B. 



r,{7r) 
22 



r,{K) 
22 



<5L5(m)=0 



+ 0{q'/N, 



C) 



1 

V2 



Ml, 



1 - 



2 5^{^Ji)Ml^ ■ 
P 



{22) 



(23) 



The second line in (|22D is easily obtained by multiplying the numerator and denominator 
with the factor (1 — 2 6Lg{fi) My^/ p). After introducing this definition, 6Lg{n) shuffles from 

i?22^'' to the parameter My^(/i). Thus instead of two independent constants, SLKfi) and My^, 
we only have the combination My_^{fi) replacing everywhere the parameter My_^. The term 

6Lg{fi) disappears from -822^^ hence leaving in the regularized Feynman integral an explicit 
dependence on /i. Moreover, the use of My^ (/x) in (|^) allows us to recover the whole value 
of the xPT running coupling 



+ 5Ll{fi) 



2Ml^{fi) 



(24) 



up to the considered order. Therefore the parameter My^ (/i) captures the right dependence of 
Lg(/i) on the renormalization scale. In our phenomenological analysis, we will adopt the usual 
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reference value /xq = 770 MeV. Later on we will perform numerical studies at different scales 
fi and will examine the corresponding values of LK/j,) derived through (plf). The prescription 
of ehminating 6Ll{jj,) from -822^'* is assumed in the following. 

When studying the experimental data, we will observe that the couplings Fy^ and Gy^ 
are not exactly the ones provided by the large- A'^c' limit, but they have small deviations. 
These parameters suffer also slight variations when more than one resonance is taken into 
account. In that case, the scale dependence does not go in such a straightforward way to 
the parameter My^(/i) as we have seen in (PBf), although the relation is still obeyed within 
a given accuracy. The other parameters are going to suffer very tiny modifications with the 
scale but, at the precision of our study, they remain like constants. 



5 Phenomenology 



We are going to analyze the experimental data for the vector form factor, which is much 
cleaner than the one from vrvr scattering. The vector form factor can be experimentally tested 
in the photoproduction of pseudoscalars from e+e" annihilation or in r decay. Although there 



are many data from e+e |20, UM, we have decided not to consider them, as we have not 



taken into account the ui-p mixing. We have studied the r Vr'^ix data from ALEPH |]TP[, 
which provides a covariance matrix to account for experimental error correlations. Similar 



data from CLEO and OPAL ^6] are also available. 



The range of validity up to which we will extend our fit is at most ^/(p < 1.2 GeV. Beyond 
this energy, multiparticle channels become important. First we perform a fit to the modulus 
of the VFF (ALEPH data) with the p(770) resonance only. This yields the parameter My^ (/x) 
and the couplings Fy^ and Gy^- We choose as matching scale /^o = 770 MeV, take the pion 
decay constant f = fn = 92.4 MeV as an input, and fit the region 2m^ < ^/s < A^^^ = 
1.2 GeV. We obtain the values shown in Table |l|, with a x^/dof= 24.8/25. The corresponding 
VFF is shown in Fig. In order to estimate the systematic errors, we have varied the chiral 
parameter / in the interval / = 92.4 ± 1.0 MeV and the final point of the fit A^^^ between 
1.0 and 1.2 GeV. All these effects yield a more conservative result with a broader error. The 
first error in Table |1| is the one provided by MINUIT while the second is our estimated 
systematic uncertainty. 

Besides the lagrangian parameter My^{fi), we can determine the more usual "physical" 
masses: the Breit-Wigner mass and the pole mass M^^". The energy where the phase- 
shift 0,r7r = f defines the Breit-Wigner mass M^^^ and the corresponding width is given by 

1/T = ^ The complex pole of the observables in the second Riemann 

^~ BW 

sheet, s^"'" = [M^"^" — irp°'''/2)^, defines the alternative pole parameters. In Table |I| we 
have written the resulting values for these two different mass and width definitions. In 
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Figure 4: VFF fits to the r —>■ Vr'^ir ALEPH data [^], with one and two p resonances. Also shown 
are the e^e^ 2tt data points from NA7 [pO| ]. 



order to derive those numbers, we have taken into account the correlations among the fitted 
parameters My-^(yUo), iVi and Gvi- Owing to the off-shell behaviour of the denominator, 
the pole mass turns out to be lower than the Breit-Wigner mass, in agreement with former 



works p8| . The opposite behaviour would have been obtained from a constant Breit-Wigner 



width parameterization. 

In Fig. ^ we plot the phase-shift 0,^71 • In the low-energy region ^/s < 0.7 GeV, the 
experimental data appears to be slightly above the predicted values. The same behaviour 
can be observed in previous theoretical studies |]14|, |16|, |1^ ^ The experimental errors 
are probably underestimated in this region, although higher-order chiral corrections could 
induce small variations to our predictions. Other studies seem to have a better control 
of the region closer to the vrvr threshold and dominated by the xPT constraints. Beyond this 
region the agreement of our one-resonance analysis with the scattering data is good up to 
y/q^ < 1 GeV. Above this point the prediction for the scattering amplitude begins to fail. 

In order to better study the region around y/s ~ 1 GeV, we include a second vector 
multiplet with the p(1450). The effect of the tail of the p(1450) can modify slightly the 
distribution in this region, where still the p(770) dominates. Nonetheless, we cannot study 
energies much higher than ^/s ~ 1.2 GeV, since some not well-known strong inelasticities do 
arise (the experimental phase-shift data does not seem to pass through 37r/2 at the p(1450) 
mass |]HT]]). Clearly, the two-pseudoscalar loops cannot incorporate all the inelasticity needed 
to describe the p(1450) region. Other multiparticle intermediate states may be responsible 
for this large effect. 

We have fitted our theoretical determination of the scattering amplitude with two res- 
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Chiral Coupling 


PI' 'Uj 


p[i (0) + p(^i4oU ) 


MvM 


845.4 ± 1.1 t^l MeV 


839.4 ± 1.41^1 MeV 




1.696 ±0.008 t^Q^oll 


1.669 ±0.008 ±0.017 


\GvJf\ 


0.695 ±0.004 t^of^l 


0.670 ±0.005 t^o^oll 


FviGvi/P 


1.178 ± 0.010 i|]:[}[J^ 


1.119 ± 0.012 


Tr( ., \ _ Y- ^ViGvi 


(7.04 ± 0.05 llii?)- 10-3 


(6.79 ± 0.09 Igi?) • 10-3 




776.0 ± 1.6 ig;? MeV 


773.9 ± 2.0 MeV 


r 

BW 


156.2 ± 1.6 MeV 


150.2 ± 2.0 MeV 




764.1 ±2. 7 t|j MeV 


770 ± 3 ± 3 MeV 


"P^polc 

p 


148.2 ± 1.9 MeV 


137.3 ±2.6 ±2.6 MeV 



Table 1: Determination of some RxT and xPT couplings, at the scale hq = 770 MeV, from the 
VFF fit. The parameters Fy-^j f and Gy^j ^ have the same sign as Fy^Gy^j p > 0. 

onances to the experimental phase-shift in the region 0.7 GeV< ^/s < 1.2 GeV. The fit is 
not very sensitive to the p(1450) mass, allowing a wide range of values. Nevertheless, it 
requires that My^^fxo) > 1550 MeV. Taking My^^po) = 1550 MeV, the fit to the phase-shift 
gives My^ifxo) = 841.8 ± 0.6 MeV, GyJ f = 0.6631 ± 0.0027 and GyJ f = 0.373 ± 0.028, 
with x^/dof= 18.8/22. The fitted value of Gy^/f increases for larger masses of the p(1450) 
resonance; the central value grows to 0.57 for My^^fi^) = 2000 MeV. The precision of Gy^/ f 
is improved, as expected, because the phase-shift has a larger sensitivity to this parame- 
ter. The differences between the analyses of ^tttt with one and two resonances are tiny for 
^/s < 1 GeV. Beyond y/s ^1.2 GeV, the description breaks down because the pathological 
^{2n ± 1) behaviour of the phase-shift in the neighbourhood of the p(1450) still remains [see 
Fig. i. 

We have performed next another fit to the VFF ALEPH data, with two vector multiplets 
and taking A^^^ = 1.2 GeV. Since in this region the data have very small sensitivity to the 
p(1450) mass and the coupling Gy^, we introduce as an input the value of My^d^o) and the 
corresponding coupling Gy^ obtained from the phase-shift fit. The results of this VFF fit, 
given in Table |I], have a x^/dof= 14.7/24. The systematic errors have been estimated varying 
the pion decay constant in the interval / = 92.4 ± 1.0 MeV and the value of My^lfio) in the 
range| from 1550 to 2000 MeV, which implies Gyjf = 0.37 ± 0.03 tg;^. In this analysis we 

^ Notice that the one-resonance results indicate that Afy(/io) is around 100 MeV larger than A/j^,,,, or 
MP°'°. The experimental situation of the p{1450) is rather unclear and it might be possible that it has an 
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Figure 5: Phase-shift (j)^T^ of the vrvr scattering amplitude. The parameters employed for the one- 
resonance graph are the same than in Fig. ^. The curve with two resonances takes the values from 
the fit to the scattering amplitude, with the inputs Mv'2(//o) = 1550 MeV and / = 92.4 MeV. 



have recovered as well the Breit-Wigner and pole masses and widths for the p(770) meson. 
We have not tried to determine the p(1450) pole, because it would lie in a region which is 
not well described. We also give in Table |l] the xPT coupling L'^^^iiq) at the matching scale 
po = 770 MeV. 

The VFF fit is sensitive to the product of couplings Fy^Gy^j p . One gets, 

Fy.GyJf = 0.007 ± 0.024 (25) 
For the range of Gy^/ f values quoted before, this implies Fy^/ f = 0.02 ± 0.06io;o8. 

Modifications of the p(1450) inputs produce sizable variations on the p(770) couplings. 
Thus, a better knowledge of the p(1450) is needed to get more accurate values of the p(770) 
parameters from a two-resonance fit. The results are consistent with the more precise deter- 
minations from the fit with only one resonance, which we take as our best estimates. 



5.1 Running of Lg{fi) 

We have seen in Section 4, from a simplified theoretical analysis, that the parameter My-^ (p) 

r ( P) 

depends on the xPT renormalization scale adopted in the loop function B22 , in such a way 

even higher mass or that a strong interference of two vectors, p(1450) and p(1700), is needed to properly 
describe the data |^ 
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Figure 6: Comparison between the usually quoted value of the chiral coupling Lgi^J-) |^ (shadowed 

band) and some determinations from the fit at several matching scales: LgifJ-) = J2i 2M^ (^) • '^^^ 

smaller error intervals are the statistical uncertainties given by MINUIT, the larger ones indicate 
the total errors including systematic contributions. 



that the physically measurable VFF is scale independent as it should. The dependence of 
(/i) with the scale was given by the equation 

M^^(/.,)-M^^(/.0 = ^ln(|), (26) 
1 1^ /"/ft 



as Ll{jj,2) — Lg(/ii) = (5^9 (/i2) —6Ll{ni) = —j^^ln yj^j- The theoretical running of Myj(/x) 
induces a scale dependence on the predicted value of Ll{fi) in Eq. (P^). When the phe- 
nomenological fit is performed at different values of fi, the parameter My^ (fi) increases with 
/J,. The other parameters of the fit remain essentially unaffected, i.e. they suffer modifications 
much smaller than their errors. Varying the scale /i in the range between 0.5 GeV and 1.2 
GeV, the varies less than 2%. 

The fitted L'g{fi) results are compared in Fig. || to the usually quoted values 0. At the 
standard reference scale /xq = 770 MeV, we obtain 

L^(/io) = (7.04±0.05l°i?)-10-3, (27) 



which improves considerably previous determinations 0, |33|. The systematic errors would 



increase to 11q[Io if we would have considered the fit with two resonances. The lack of knowl- 
edge about the second multiplet parameters introduces an extra uncertainty of the same 
order than the one we have with only one resonance. 
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Figure 7: Allowed region in the {Rv,Ra) plane. The shaded zone on the left corresponds to the 
one-resonance study, while the dotted zone results from the analysis with two resonances. 

5.2 Ijar^e—Nc relations 

As we work at higher orders in I/Nq, our experimental results have next-to- leading devia- 
tions from the LO values provided by the two short- distance QCD relations (|^) and (|^). We 
are going to test now how well they are satisfied. Typically, there should be a deviation from 
zero of 0{1/Nc) in the VFF relation (with Nc = 3 in physical QCD), as the leading terms 
of the left-hand side of the equalities are Fy^Gvi/ P ~ 1- The deviation in the axial form 
factor constraint should be of 0{1/Nc) ■ 0.03, because its leading terms are 2Fv^Gv^ I My^ ~ 

In Fig. ^, we have plotted the variables 

(28) 

R'^ - {e ^''''•%^- } ■ (0.03)-' . 

which have been normalized with appropriate factors so that the expected deviations from 
zero are of 0{1). We have performed a scanning of the range of values for the RxT couplings 
obtained from the VFF fits. We can see in the figure that the separation from the large- 
Nc QCD relations is indeed of the expected order for both types of fits (with one or two 
resonances). Thus, the short-distance relations (0) and @) are well satisfied, within the given 
accuracy. 
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6 Uncertainties from Higher-Order Corrections 



There exist many more diagrammatic contributions which have not been included in our 
results. We show in Appendix C that, when the production of multiparticle states is ne- 
glected, it is possible to define a generalized summation of Feynman diagrams with two-body 
topologies. It makes use of a kernel function /C, associated with the two-body scattering 
amplitude, which incorporates those contributions not included in our effective s-channel 
vertex of Fig. H(b). The resulting VFF can be formally written in a very compact form, 
given in Eq. ( |C.10| ). Making a 1/Nc expansion of the kernel /C, one can easily check that 
our s-channel result in Eq. (|l^ corresponds to the leading-order approximation. The first 
correction originates from a single resonance exchange in the t-channel, which induces a sub- 
leading contribution of 0{1/Nc) to the kernel /C. The exchange of n meson fields contributes 
to the kernel at 0(1/Arg). 

A general calculation of those higher-order corrections is a formidable task. We know, 
however, that in the energy region we are studying the tree-level scattering in the t-channel 
is much smaller than the one coming from the s-channel, what seems to imply that they 
contribute as a small perturbation. To estimate the size of those corrections, we have ana- 
lyzed the leading contribution from t-channel resonance exchange between the final pions. 
According to the results in Appendix C, it induces a multiplicative correction into the VFF: 

^0, (29) 

where is the contribution from a single t-channel exchange. 

The complete calculation of is rather involved, since it makes necessary to address 
the renormalization of R^T. This is a very interesting issue, which we plan to analyze in a 
future publication where a full analysis of the VFF at next-to-leading order in 1/Nc will be 
attempted. Here, we are only interested in its numerical impact on the results presented in 
the previous sections. For simplicity, we will study Q\j^ in the SU(2) theory; i.e. we neglect 
the tiny contributions from diagrams with kaons in the intermediate loop or in the final state 
^gt^^,K) _ gt^^K,TT) _ gt^^K,K) _ Moreover, we will work in the chiral limit (m^ = 0). 

Although there are several Feynman diagrams contributing, we only need to consider 
the dominant one where the current vertex {ki — k!^) produces a vr~(A;i)7r°(A;2) pair, which 
is rescattered through a t-channel resonance. This diagram generates the interesting non- 
analytic contributions, plus a divergent local correction which should combine with the local 
contributions from the other diagrams to provide a physical finite result. Since we are in- 
terested in the region ^/s -C 2Mv (i.e. we work below the two-resonance cut), there are 
no additional sources of non-analytic terms. The local ambiguity can be fixed to 0{E^) by 
matching Eq. (|29|) with the known xPT result. This requires Q\j^ ~ (9(i?^). The exchange 
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of a vector resonance can be easily computed in this way. One gets: 



Li2 1 + 



+ In 




- Li2 (1) 



+ 4 + 



6 M, 



r- Ml 



9 
2 



35 



36 M, 



(30) 



Since g\ 



we can neglect the exchange of higher- mass vector resonances. 



However, we will also consider the t-channel exchange of scalar resonances from the lightest 
multiplet, with a mass Ms — 1 GeV ||10|, ^ and couplings Q, — //2 |^. It provides the 
contribution: 



y^t (7r,7r) 



Lio 



+ ln 



24+2 

. 1 1 . 



(31) 



Mi 



6 



1 1 
~ 2 ^ 36 



This result includes contributions from the singlet and the octet scalars. 

At energies below and around the p(770) peak, these t-channel diagrams give a correction 
smaller than 5%, which is within the uncertainties of the numerical analyses performed in 
the previous section. However, above ~ 1.2 GeV the vector contribution becomes larger 
than 10% and these topologies cannot be neglected any more. This kind of diagrams turn 
out to be very important at high energies. 



We have repeated our previous fits to the VFF ALEPH data, including the correction 
induced by Q\^'^'^ = Q\v''^'^ + Q\s'''^\ The results of these fits are compatible with the ones 
obtained before, showing that our former studies neglecting crossed channels provide a good 
description within the given precision. 



7 Conclusions 



A quantum field theory description of strong interactions at energies around the hadroniza- 
tion scale, £" ~ 1 GeV, requires appropriate non-perturbative tools. While a fundamental 
understanding of the confinement region of QCD is still lacking, substantial phenomeno- 
logical progress can be achieved through effective field theories incorporating the relevant 
symmetries and dynamical degrees of freedom. 

Using an effective chiral lagrangian which includes pseudoscalars and explicit resonance 
fields, we have investigated the VFF and related I = J = 1 observables in the interesting 
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~ 1 GeV energy range. The heavy particles make the standard chiral counting in powers 
of momenta useless, because their masses are of the same order than the chiral symmetry 
breaking scale. Therefore, we have adopted instead the more convenient large- A^c- expansion, 
which provides a powerful tool to organize the calculation. 

At the leading order in 1/Nc, one gets an excellent description of the VFF, far away 
from the resonance singularities. A proper understanding of the zone close to the p(770) 
pole, requires the inclusion of next-to-leading contributions providing the non-zero width of 
the unstable meson. The dressed propagator can be calculated through a Dyson-Schwinger 
summation of the dominant s-channel rescattering corrections, constructed from effective 
Goldstone vertices containing both the local xPT interaction and the resonance-exchange 
contributions 



13, 14, 15, 16 



We have extended the Dyson-Schwinger summation of effective vertices to handle prob- 
lems with coupled channels in a systematic way, through the recurrence matrix Ai. The 



generated by final-state interactions, provides the right unitar- 
0|. Moreover, with an S'?7(3)-symmetric dynamics 
'1 — A1)~^ acts just like a pure 



inverse matrix {1 — Ad)' 
ity structure of the observables |]l^ ^ 
(the vertices contain only derivatives and no quark masses 

number [1 — tr{7Vl}]~^. Hence, there is no mixing among loops and the total decay width 
is simply given by a sum of separate contributions from the different channels, which corre- 
spond to the partial decay widths. An improved diagrammatic summation of more general 
two-body topologies has been given in Appendix C. It includes the smaller t-channel correc- 
tions, through the 1/A^c expansion of a non-trivial interaction kernel JC associated with the 
two-pseudoscalar scattering amplitude. 

The Feynman loops fully determine the non-analytic contributions, which are dictated 
by unitarity and chiral symmetry. The local corrections, however, are functions of the theo- 
retically unknown couplings of the effective lagrangian. They incorporate the short- distance 
dynamics and take care of the regularization and renormalization prescriptions adopted in the 
calculation. A significative reduction on the number of free parameters is obtained, requir- 
ing the different amplitudes to satisfy the appropriate QCD constraints at large momentum 
transfer [^, In fact, a very successful prediction of the most relevant 0{E^) xPT cou- 
plings is obtained, under the reasonable assumption that the lightest resonance multiplets 
give the dominant effects at low energies |]13[. We have resolved the local ambiguities of the 
VFF, imposing the QCD short-distance constraints and performing a low-energy matching 
with the known 0{E'^) xPT result. 



T3|, we have obtained a good fit 



in the range 2m^ < V¥ < 1-2 GeV. At the chiral 



Working within the single-resonance approximation 
to the ALEPH r z/^27r data |g, 
renormalization scale /xq = 770 MeV, the fit gives the values shown in Table |l| for the main 
p parameters. The corresponding resonance pole 5^°''' = (M^"'" — irp°'''/2)^ in the second 
Riemann sheet is found to be at: 

ppoic ^ 2 ± 19 +1.7 ^ (32) 



M;°'" = 764.1 ± 2.7 tti MeV , 
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We have achieved an improved determination of the xPT couphng 



L^(/io) = (7.04 ± 0.05 ^[Ji?)■10-^ (33) 

at fit) = 770 MeV. Performing the phenomenological fit at several scales fi, ones obtains the 
proper running of L^fi) as prescribed by xPT. 

To test the convergence of the 1/Nc expansion, we have analyzed the deviations between 
the fitted parameters and the corresponding theoretical large- A^^c predictions [Tl|. The dif- 



ferences are found to be of the expected 0{1/Nc) size, showing that the limit A^^^ oo 
provides indeed an excellent description of the local chiral couplings. 

We have also investigated the corrections induced by the tail of the p(1450) vector res- 
onance at the higher side of our energy range. The effects are sizable, but the sensitivity 
is not good enough to make a precise determination of its parameters or to disentangle the 
existence of several higher-mass states. In order to do that, one would need to study higher 
energies where other multiparticle final states, beyond the two-body modes that we have 
analyzed, become relevant. Moreover, a better calculation of t-channel contributions would 
be needed, because they are no longer small above 1.2 GeV. 

To summarize, we have performed a detailed analysis of the p(770) region, imposing all 
known theoretical constraints. The main p parameters and the xPT coupling LKfi) have 
been determined with rather good precision. More work is needed to extend the results at 
higher energies. It would also be very interesting to investigate in a similar way the scalar 
sector, specially the pathological I = J = observables. We plan to address these issues in 
forthcoming works. 
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Appendix A: Feynman Integrals 

(p) 

The loop function used in the text is defined through 



J i(27r 



dk'^ 



{2t:Y {k^ - ml) [{q - kf - ml] 



with 



B. 



(P) 
22 



1927r' 



1 - 



6m I 



Aoo + In 



ml 



= -D22 Q 9 + -D21 Q Q , 



8mp 5 o , / ap + 1 

+ o+4ln - 

q 3 \ap — 1 



(A.l) 



(A.2) 



where Ago = 7^ /u'^ + 7^: — In (47r) — 1, je — 0.5772, is the renormahzation scale and 



ap = — 4mp / is the usual phase-space factor. 



The real part of this Feynman integral is divergent, but its imaginary parts is finite and 
takes the value 



1927r 



(A.3) 



The dilogarithm function which arises in the crossed-channel calculations is defined as 



Li^ 



iiy) = - I — ln(l - xy) = - / — In(l-x). 
J X J X 

It has an imaginary part given by 

Im{Li2(y + ie)} = vrln (y) ^(y - 1) . 



(A.4) 



(A.5) 



Appendix B: Matrix Relations 

In the isospin limit, the matrix E~^T^q is proportional to a dimension-one projec- 
tor. Therefore, it obeys the properties of a general dimension-one projector P and a general 
matrix B: 

P-B-P = AP, (B.l) 

with A = tr{P ■ B}. When the inverse matrix (1 — P • B)~^ is multiplied by the eigenvector 
u of P, or by the matrix P, we obtain 



(l-P-B)-^ it 



1 



u . 



{1-P-B)~^ P 



1-A ' 
' P. 



1 - A 



(B.2) 
(B.3) 



In the study carried on before in Section 3, the matrices P and B were S T^q E and 
(— 1927rP22), respectively. The matrix P • B was just A4 and the vector u was JFq. 
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Appendix C: Summation of General Two-Body Topologies 



The Dyson- Schwinger summation performed in Section 3 incorporates the dominant 
s-channel contributions. Moreover, the adopted matching procedure to the low-energy xPT 
results takes care of tadpoles and local contributions, to the considered order in the mo- 
mentum expansion. There are, however, many more diagrammatic topologies which have 
not been considered yet. Neglecting the small corrections coming from multiparticle inter- 
mediate states, it is possible to define a generalized summation of Feynman diagrams with 
two-body topologies. 

As we saw before, the effective vertex in Fig. ^(a) for the vector current insertion producing 
a P~P^ pair of pseudoscalars shows the momentum structure: 



V2, 



(C.l) 



with P^:'^ 



.AlV 



^ and P^ 



the usual transverse and longitudinal Lorentz projec- 
tors. In a similar way, the effective vertex in Fig. H(b) describing the s-channel scattering of 
two pseudoscalars, when projected on the P-wave (/ = J = 1), takes the form: 



487ri 



{h-k2U, (C.2) 



with Pi, P2 {ki, ^2) the outgoing (incoming) momenta. The matrix T£*q is the corresponding 
I = J = 1 partial-wave scattering amplitude. 

Let us define a general kernel /C*^'^'"^(/ci, ^2,^1,^2) associated with the two-body scattering 
amplitude from (n)-type pseudoscalars to (m)-type pseudoscalars. This kernel, shown in 
Fig. 0(c), contains the identity operator (no scattering) plus all interaction diagrams without 
intermediate effective vertices 



Now let us connect the effective vector current insertion Vf^ to the kernel /C*^™''") , as shown 
in Fig. 0(a). The outgoing pseudoscalars from the kernel are joined again into an effective 
scattering vertex ( |C.2| ). This generates the dressed structure: 



{Pi - P2Y 



^-^T'l^{s)^~^ Pl% 



[n(s) PtI + n'(s) PlI] 



fo{s)PT^, + f'o{s)PL^^]V2, 
with the matrices n(s) and n'(s) defined from the kernel integral 

n("^.")(s)PrX + n'('™'")(s)Pr:: = f-^i X 



(C.3) 



T^tF^ (^'1 - k'2)aA^'-\k'l)A("^\k'l)IC^"^'-\k\, k'2, fci, k2)A(-\kl)A(-\kl){k, - k,r 
(27r) (27r) 

(C.4) 
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Figure A: Basic pieces of the general summation of two-body topologies. The first row shows the 
general kernel /C while the second one only includes the contributions from ladder diagrams. 



where A*^™^(A;^) is the propagator of a (m)-type pseudoscalar. Performing the trivial products 
of Lorentz projectors, the vector current matrix element with one intermediate kernel and 
ending into an effective scattering vertex takes the form: 

= {p^-p2y[T.-'Tl^{s)T.-'li{s)fio{s)PT^, + E-V^Js)E-^n'(s)A(s)PL(;] . 

(C.5) 

When the outgoing pseudoscalars are both on the mass shell the longitudinal term becomes 
zero. 

We can easily iterate this algebraic procedure and consider a series of N intermediate 
kernels and effective scattering vertices, attached to the current insertion. The first kernel is 
connected directly to Vq] then it comes an effective scattering vertex %, followed by another 
kernel, and so on. The outgoing pseudoscalars are attached to the final effective vertex. The 
resulting contribution to the VFF is expressed as: 

vj^ = (pi -P2Y [ (s-^T^^oE-in)^ Pr'i + (s-^T'^oE-^n')^ A PL'i]V2. (C.6) 

Thus, the summation from = to infinity becomes: 

V>' = (Pi -P2y [(1 - E-iT,^oS-'n)-i Prt + (1 - E-^T'^oS-^n')-^ /'o PLt]V2. 

(C.7) 

This sums all diagrams ending in an effective scattering vertex. Finally, we add the 

diagrams where the last effective vertex is connected to the outgoing pseudoscalars through 
the kernel. This extra contribution is given by the form factor of the factorized element 



^2,Pi,P2) A(-)(^?) lsS^\kl) {h - k^r , 



(C. 
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shown in Fig. 0(a), which we have separated into transverse and longitudinal parts. The 
summation of all types of diagrams gives then, 

= {P1-P2Y [Q{1 - T.-^Tlj:-^Tl)-^ Pt^ + Q'{1- T.-'T'lj:-'Tl')-' /'o Pl':]V2. 

(C.9) 

With the outgoing pseudoscalars being on-shell, the resulting VFF takes the compact 
form: 

= g ■ (i -^-^Ti^^'^ny^ ■ pQ. (c.io) 

The simplest kernel is the trivial direct connection of the incoming and outgoing pseu- 
doscalars (/C = /). In that case, the integral ( |C.4| ) reduces to the usual two- propagator loop. 



n = —192tc B22, and Q = I. One recovers then the expression (14), obtained through a 
Dyson-Schwinger summation of s-channel scattering vertices. Eq. ( |C10D provides a system- 
atic way of improving the result, with the use of more complex kernels. The calculation could 
be organized with the use of a 1/Nc expansion of the kernel /C; the trivial identity opera- 
tor corresponding to the lowest-order approximation in this expansion. The first correction 
comes from a single resonance-exchange in the t channel, which induces a contribution of 
0{1/Nc) to the kernel. The exchange of n meson fields would contribute at 0{1/Nq). 



C.l Ladder Diagrams 

The calculation of higher-order diagrams with an arbitrary number of resonances exchanged 
in the t channel turns out to be a very complicated problem as each loop is connected to 
others. However the optical theorem relates the form factor diagrams Fig. 0(d) with the 
scattering amplitude through ladder diagrams. Fig. 0(f), in the familiar way [jT^ ^ 



implying 



Im T* 

im 



. ^3 . gt* 



(C.ll) 
(C.12) 



rpt 



T*-^ + 0(l)-^S, 

i + o{i/Nc)-i{j:-X 



-Irr^t + C(l/iV^)l 



(C.13) 
(C.14) 



where the terms C(l), 0{1/Nc) and 0{1/Nq) correspond to NLO contributions in l/Nc, all 
of them real in the physical region when multiparticle channels are neglected. The matrix T^q 
is the tree-level scattering amplitude through a crossed resonance and the diagonal matrix 
is just the phase-space matrix but with each ap multiplied by a threshold factor 9 {q"^ — Ami,). 

The basic behaviour of these quantities is driven by the tree-level term, as the crossed 
scattering amplitude is tiny at the energies we are considering. It turns to be important at 
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very high energy, where the t-channel becomes the dominant amphtude. Thus, the matching 
of to the lowest-order contribution plus the diagrams with only one t-channel resonance 
exchange is a suitable assumption: 

(C.15) 

with Im = S-^T* ■ S3. 
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